Structural instability of vortices in Bose-Einstein condensates 
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In this paper we study a gaseous Bose-Einstein condensate (BEG) and show that: (i) A minimum 
value of the interaction is needed for the existence of stable persistent currents, (ii) Vorticity is not 
a fundamental invariant of the system, as there exists a conservative mechanism which can destroy 
a vortex and change its sign, (iii) This mechanism is suppressed by strong interactions. 
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Superfluidity is an intriging concept which dates back 
to the discovery of He-II Q , and which is by now text- 
book material. Among the properties of a superfluid 
there are irrotational flows, quantization of circulation 
and stability of persistent currents. The theoretical dis- 
advantage of He-II is that the strong interactions which 
take place in the liquid prevent a simple microscopical de- 
scription of the system, and make it impossible to clearly 
untangle the ultimate cause of superfluidity. 

The achievement of Bose-Einstein condensation with 
ultracold atomic gases has revived the interest on a the- 
oretical description of superfluidity. The flrst reason is 
that Bose-Einstein condensation has long been thought 
to be the cause for the superfluid behavior in He-II [||. 
But the second reason is that gaseous condensates can be 
accurately described with simple mean fleld models that 
can be subject of both analytic and numerical study |3|. 

In the last few years, the superfluid nature of these con- 
densates has got experimental demonstrations by several 
means. One of the most striking features of superflu- 
ids is the quantization of circulation, which leads to the 
existence and stability of vortices even with the lack of ex- 
ternal forcing. Vortices have been theoretically predicted 
to be present in gaseous condensates with repulsive in- 
teractions |^-||], a prediction which has been conflrmed 
in recent experiments ]9|-pj]] where vortices have been 
produced and shown to be, at least, long-lived. 

The stability of vortices as seen in experiments and in 
numerical studies, has lead us to consider vortices and 
their topological charge as robust entities. In this work 
we prove the surprising fact that a weakly interacting 
condensate in a stationary trap cannot host stable per- 
sistent currents of well deflned vorticity, showing that 
nonlinearity is therefore essential to sustain vortices. We 
show what is the mechanism for the structural desta- 
bilization of vortices and propose a way to verify this 
prediction with in current experiments. 

The model.- In this work we study a dilute gaseous 
BEG in the zero-temperature limit. This system is ruled 
by the Gross-Pitaevskii mean fleld equation H 
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where A = X]j=i^^/^^j n-dimensional Lapla- 

cian, V(r) = ^muo'^ X]j=i external potential 

which conflnes the condensate and is 5-wave scatter- 
ing length for the binary collisions within the condensate. 
We will focus on condensates with > 0, which repre- 
sents the repulsive interaction between bosons. The norm 
of the wavefunction corresponds to the number of bosons 
in the trap, N = 

Due to the quantum nature of the system we are study- 
ing, the velocity of the fluid can be derived from the 
macroscopic wavefunction of the condensate 
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Here we have separated the wave function into its mod- 
ulus, 1^1 = y^, and its phase, = arg^. 

From now on we will use adimensional quantities to 
ease the analysis. To do so we deflne a new set of units 
based on the trap characteristic length, ao = ^/h/m^J^ 
and period, T = 1/uj deflned as Xj = rj/ao^t = r/T and 
a new wavefunction '^(r, t) = ^(x, r) satisfying 
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where U = Alias /clq- The value of NU measures the rel- 
ative importance of the nonlinear term /7|?/^p. Through 
this paper we will call a weakly interacting condensate a 
system in which NU < 1 as oposed to a strongly inter- 
acting condensate in which NU ^ 1 [|l2| . 
The energy of the system is given by 
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which is a conserved quantity. 

In this work we will consider a BEG placed in a tighly 
conflning pancake type trap [|l3[, for which the proflle 
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of the solution along the confined dimension is very well 
approximated by a gaussian function [Q . The procedure 
for reducing the Gross-Pitaevskii equation to an effective 
two-dimensional equation is described in detail in Ref. 
]T5[ . The only relevant difference with Eq. (||) reduced 
to two spatial dimensions is that the nonlinear term is 
reduced by a factor related to the trap asymmetry. We 
will consider then that our effectively two-dimensional 
system is described by Eq. (||) but with a new U = 
aAnaso^/ao [ p!6[ . 

Noninteracting systems.- A persistent current is a 
stationary solution in which the velocity field is nonzero, 
V 7^ 0, and remains still, 9tv = 0, even without external 
forcing. Therefore, any solution of Eq. (|l|) corresponding 
to a persistent current must have a position dependent 
phase. To ease the analysis we will write such stationary 
solution as ?/^^(x, t) = 0^(x)e*^^ = + ^'^m,*^ where 
tjj^^r and ip^^i are real functions. 

We will first consider the noninteracting case {U = 0), 
in which the persistent current satisfies 
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with the Hamiltonian H = — + ^(x). Remarkably, 
Eq. (H) implies that ijj^^r and ip^^i must be eigenstates of 
H with the same eigenvalue. Therefore, in noninteracting 
systems, persistent stationary currents can exist only if 
the spectrum of the Hamiltonian is degenerate. 

Furthermore, degeneracy is a very specific situation 
which can be broken by almost any perturbation -e.g. 
any small asymmetry of the confining potential-, which 
leads us to conclude that persistent currents in nonin- 
teracting systems are structurally unstable^ a term which 
will be precised later. 

To understand what happens to the quantum fiuid 
in the noninteracting regime, let us consider a simple, 
two dimensional harmonic oscillator whose three low- 
est energy eigenstates of H are exp (^—f^ ~ ^) ' 
(j)x oc X(/)o, (l)y cx ?/^0 7 where we make use of the adi- 
mensional natural widths of the oscillator, a = 

b = i/^. 

In the case without interaction {U = 0) and ra- 
dial symmetry (a = 6), vortices are solutions of 
the type tlj{x^y^t) = \/N/2 [^^^ + i^^/] ^^^^^ which 
is an eigenfunction of H. The analogous of 
this current for the asymmetric case (a ^ b) is 
iIj{x, y, t) {x^/N^ + iy^/Jlye^'^^y~^'^^^) e*^^^0o- Due to the 
existence of a time dependent phase between both com- 
ponents, the solution is not stationary and evolves far 
from the initial configuration. In fact, for t = tt/ {iiy — fix) 
the relative phase leads to a complete inversion of the 
topological charge of the vortex from m = +ltom = — 1 
passing through an intermediate state where a dark line 
breaks the condensate into two pieces. 

Structural stability.- Thus, for the noninteracting 
case: (i) There are no stationary vortex solutions in any 



stationary asymmetric trap (all eigenfunctions are real), 
(ii) Vortex type solutions are unstable in those traps since 
they follow an evolution which leads to an inversion of the 
topological charge. 

The vortex reversal can be induced by any small per- 
turbation that breaks the degeneracy, such as small 
asymmetries or even spatial noise. This is an example of 
structural instability. The important difference between 
stability referred to perturbations of the initial conditions 
and robustness referred to perturbations of the governing 
evolution equations, is not always properly appreciated, 
but it is a key ingredient of the phenomenon discovered 
here. 

It is important to stress that our previous finding im- 
plies that there exists a fundamental limitation on the 
stability of vortices, which makes that as soon as the 
forcing is removed, these currents will persist only for a 
period of time that depends on the strength of the inter- 
actions and the amount of perturbation. 

Role of interactions.- Let us now consider the inter- 
acting case and look for the simplest stationary currents 
on a two-dimensional condensate in the very weakly in- 
teracting limit. In this limit the ground state of the con- 
densate is almost equal to the ground state of the linear 
problem ?/;o, and for the vortices we may use the simple 
ansatz = ^/N^cj)^ + i^fNy(j)y^ where and Ny are 
variational parameters and can be regarded as the num- 
ber of atoms on each dipole state. The energy of our 
ansatz is obtained substituting into Eq. (H) 

M.iV. + ^iyNy + ^ (3iV| 
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where jj^x = fA^^ + ^A^y and /i^ 
eigenvalues of the (j)x^y basis functions. 
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FIG. 1. (a) One of the dipole solutions in a very weakly 
interacting condensate, (b-c) Density plots of stationary vor- 
tices for UN ~ 14, 31 respectively. All plots expand the spa- 
tial region [—3, 3] x [—3, 3] and correspond to a trap with 

We have to search the critical points of (||) over the 
line N = Nx + Ny. When U = the energy functional 
has only two such critical points which correspond to the 
(j)x and (j)y stationary states. When U we still find 
those two dipole solutions with energies given by 

Aiy + ^)iV, N = Ny,N, = 0. (8) 
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Furthermore, if and only if the interaction exceeds a 
critical value 



UN > 27rab\jiix — fJ^yl = 
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there exists a third solution, which is a deformed vortex^ 
and whose energy and populations are 
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According to these equations we classify the (A/", U) 
space into several regions. For UN < a, no vortex is 
found. Above UN > a not only the vortex exists but 
it is also energetically more favorable than the dipoles 
therefore preventing the possibility that the cloud spon- 
taneously splits into two pieces. 

Numerical simulations.- To study the sensitivity of 
the system to finite (but may be "small" ) perturbations 
we have proceeded as follows. First we have found the 
stationary vortex profile for given asymmetries of the 
trap and interaction values [See e.g. Fig. 0(b-c)] us- 
ing the numerical method developed in Ref. |l7|] . Next 
we have used this stationary solution as initial data for 
a numerical simulation of the dynamics of Eq. (^) in a 
perturbed system where the perturbation used consisted 
on a change of the trap parameters. 

For values of UN < a we were not able to find any 
vortex solutions, in agreement with the prediction from 
the variational analysis. 
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FIG. 2. Rupture of a symmetric vortex (Ax 



1, 



UN ^ 2.2) placed in a trap with A^c = 1, Ay = 0.77. Shownare 
density plots and phase (interference) plots on the spatial 
region [—2.5,2.5] x [—2.5,2.5] (adimensional units). 

However, we find that for medium or small UN^ where 
vortices already exist and are suposed to be linearly sta- 
ble, the application of finite perturbations leads to the 
destabilization of the vortex and the generation of a dy- 
namics with some analogues with the vortex inversion 
mechanism discused above in the framework of the non- 
interacting system (see Fig. The smaller the value 



of UN, the smaller the amount of perturbation which is 
sufficient to destabilize the system. For moderately large 
values of the self interaction structural stability adds to 
the dynamical stability that was previously known to ex- 
ist [§,||J|]. In this case vortices are stable and keep their 
topological charge even when oscillations are induced to 
the vortex profile (See e. g. Fig. ||) by strong pertur- 
bations of the confinement. In this region, neither small 
noise nor changes in the confinement are able to break 
the vortex. 




FIG. 3. Stability of a symmetric vortex (A^ = A^ = 1, 
NU ^ 30) placed in a trap with A^ = 1, A^ = 0.71 (a ^ 2.2). 
Shownare density plots and phase (interference) plots on the 
spatial region [—3,3] x [—3,3] (adimensional units). 

We have also verified that the phenomenon of vor- 
tex charge fiip is not exclusive of vortices in asymmet- 
ric traps. In fact, as the general analysis of the linear 
problem suggests, any symmetry breaking perturbation 
should lead to a simmilar phenomenon. To test this idea 
we have studied numerically the effect of a thin and in- 
tense pinning potential placed asymmetrically with re- 
spect to the trap center and verified that the same mech- 
anisms involving topological charge inversion are present. 
This interesting result adds value to the prediction since 
the study of the effect of pinning potentials on vortices is 
one of the goals of several experimental groups (see e.g. 
]rT[). In fact many other dynamical features are found 
and will be reported in detail in a future paper [Q . 

Discussion and conclusions.- In this work we have 
studied the role of interactions on the properties of vor- 
tices of a gaseous, trapped Bose-Einstein condensate. 
Our main conclusion is that vortices in a weakly inter- 
acting condensate are structurally unstable and a critical 
value of the interaction is required for the condensate to 
host stable stationary currents. 

Although the variational approach states that for 
any NU > a vortices exist and become energetically 
favourable against dipole states, our numerical simula- 
tions suggest that there exists a finite range of weak in- 
teraction values in which the vortex is close to structural 
instability and breaks under small perturbations. 

It is important to stress that this phenomenon may 
be experimentally observed using either asymmetric pan- 
cake traps |l3| or Feschbach resonances [|l9|,|o| . 

As we have discussed above, the effective interaction 
in a quasi two-dimensional system is reduced by a fac- 
tor depending on the asymmetry of the trap Using 
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the numbers of the pancake trap recently depeloped by 
MIT group to test low dimensionality effects [|l^ we get 
that U2D — 3 X 10~^. This estimate implies that a num- 
ber of particles about :^ 10^ would lead to a value of 
UN = 30 which is in the parameter range in which the 
phenomenon should be observable. It is interesting that 
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the trap used in |13| is already asymmetric since cOxf^TT = 
10 Hz and ujyl2i\ = 30 Hz, i.e a ~ 10. These numbers 
mean that the vortex topological charge inversion should 
be observable by just placing it into one such trap (using 
e.g. a rotating pinning such as the ones used by the same 
group and eliminating it once the vortex is generated). 
We have verified numerically that the vortex charge fiip 
is present in this parameter range and even for larger val- 
ues of N. The fact that in the experimental setup the 
trasverse asymmetry is not "small" helps in making the 
phenomenon accesible. 

Other possibility for observing the phenomenon could 
be to reduce the effective interacion using Feschbasch res- 
onances |l^ until a level at which the vortex would be- 
come "unstable" . This procedure may be experimentally 
controlled with high precision and very low loss of atoms 
(e.g using Rb^^ the JILA group has been able to reach 
a regime of zero effective interactions [Q). In the later 
case the preparation of a vortex should not be a problem 
since it could be generated by using either a rotating trap 
or a multiple condensate system and subsequent evapo- 
ration of one of the species ]2l|] . 

The lack of structural stability for vortices is a surpris- 
ing result. It has been widely accepted that the existence 
of a macroscopic wave-function is a sufficient condition 
for the existence and stability of vortices no matter what 
the actual model is. The dynamical stability of vortices 
as seen in experiments and in numerical studies with sym- 
metric [HJsl and asymmetric [|[ vortices in strongly in- 
teracting condensates, and some results with symmetric 
vortices in weakly interacting condensates [||, has lead 
us to consider vortices and their topological charge as 
robust entities. 

Nevertheless, as we have shown in this work there ex- 
ists a mechanism which changes the topological charge 
of the vortex while preserving the energy of the cloud. 
This mechanism is completely different from the dissi- 
pative mechanisms usually considered in BEG Q, which 
imply loss of energy and induce a spiraling of the vortex 
out of the condensate. In fact, the mechanism presented 
here is reminiscent of a phenomenon recently discovered 
in light beams , where vortices where observed to un- 
dergo extremely sharp Berry trajectories that, in the ex- 
perimental plane of observation, appear as an edge-line 
dislocation which eventually yields the inversion of the 
vortex charge. 

These predictions, striking as they are, shed light on 
the essential role played by nonlinearity in these systems 
and it is another manifestation of the richness of the non- 
linear quantum phenomena appearing in BEG. 
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